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Abstract 


The concept of local dynamical equilibrium In the form of 'eoatrophlc or 
cyclostrophlc equilibria, «^lch Is valid at high altitudes abov*' the solid' 
surface, may be adequate to describe the wide range of superrot tt Ion rates 
olfserved In planetary atmospheres. It does not address the pro )lem of 
"exactly how this superrotatlon ... Is produced and maintained gainst the 
tendency for friction to oppose differential motions between the atmosphere 
and the underlying planet" (Hide et al. , Nature. 1980). Our aialysls, to 
some extent pedagogical. Is based on the concept of a thermally driven zonally 
symmetric circulation (e.g. Leovy, 1973; Glecasch, 1975). On p'anets, with 
rotation axis perpendicular to the orbital plane, the equatorla region Is 
heated preferentially. The atmosphere expands, and large scale motions In the 
form of Hadley c«ills develop which produce a redistribution of tass and 
angular momentum. To provide a basic understanding of the phenomenon simple 
analogies are proposed, where mechanical devices with variable moment of 
inertia are frlctlonally coupled to a rotating turntable. Using a three 
dimensional zonally symmetric spectral model and Laplace transformation, the 
time evolution of a fluid leading from ^torotation under globally uniform 
heating to superrotatlon under globally ncnunlform heating Is simulated. For 
high viscosities atmospheric superrotatlon can be understood In analogy with a 
pirouette. During spin up angular momentum Is transferred to the planet. For 
reasonably low viscosities, however, this process Is reversed. The Increased 
tendency toward geostrophy, combined with the Increase of surface pressure 
toward the poles (due to meridional mass transport). Induces the atmosphere to 
subrotate temporarily at lower altitudes. The resulting viscous shear near 
the surface thus permits angular momentum to flow from the planet Into the 
atmosphere where It propagates upwards and, combined with the change In moment 
of Inertia, produces large superrotatlon rates at higher altitudes. 
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INTRODUCTION 


The phenomenon of atmospheric superrotation is ubiquitous in our solar 
system> For Earth, there exists evidence suggesting that the rotation rate of 
the thermosphere Is as much as a factor of T « 1.4 larger than that of the 
surface (e.g., King Hele, 1964). Superrotation Is observed on Mars (e.g. , 
Conmh et al. , 1973), Jupiter (e.g., Smith et al. , 1979) and Saturn (e.g. 
Smith et al., 19G1), the most striking example being that of Venus (e.g., 
Schubert et al. , 1980) where the rotation period of the visible cloud top is 
r ■ 60 times shorter than that of the solid body (243 days). Each of these 
atmospheres has a distinctive rigid shell (global average) component of super- 
rotation, Irrespective whether the planet rotates in the prograde (e.g.. 

Earth) or retrograde (e.g., Venus) direction. 

The difficulty of understanding superrotation has recently been 
emphasized in a paper by Hide et al. , 1980 from which we quote; "The density 
of the atmosphere decreases rapidly with height and more than 95% of the mass 
of the atmosphere is contained within the troposphere, which extends from the 
ground to an altitude of ~ 10 km, and the lower stratosphere. These regions 
of the atmosphere rotate faster on average than the underlying solid Earth 
with a typical relative speed of about 10 ra s”^. Exactly how this super- 
rotation (a term used for the generally more rapid relative motions found in 
the very tenous regions at heights around 300 km) is produced and maintained 
against the tendency for friction to oppose differential motions between the 
atmosphere and the underlying planet is a central problem of dynamical 
meteorology." 

Numerous theoretical models have been proposed to describe the effect 
(e.g. Schubert and Whitehead, 1969; Thompson, 1970; Gierasch, 1970; Hide, 

1966, 1971; Gold and Soter, 1971; Gierasch and Stone, 1968; Gierasch, 1975; 
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Leovy, 1973; the review of Rlehbeth, 1972) • but none 1* fully ac::epted. Vlhlle 
different processes are probably Important in different regions of the 
atmospheres and in various planetary settings, one feature all sliare in common 
is differential heating by solar radation which has long been kno»m to produce 
pressure variations and zonal wind velocities (e.g. Lorenz, 1967 )• It is the 
purpose of this paper to examine the angular momentum budget of this mechanism 
and to provide a mechanistic interpretation that nay contribute to a basic 
understanding of the average or rigid shell component of superrotation* The 
i'jiportant problem of differential rotation and equatorial superrotation will 
be discussed in a subsequent paper (Mayr and Harris, J.981). 

Leovy (1973), Gierasch (1975) and others (e.g. Tlirlkov and Challkov 
(1971), Kalvay de Rivas (1975)) proposed models for the zonally symmetric 
circulation with superrotation resulting from differential heating. It was 
suggested by Leovy and shown by Gierasch that superrotation is maintained in a 
balance between upward transport through the meridional cell and downward 
transport through eddy processes. With a steady state solution of the 
nonlinear energy, continuity and momentum equations, Gierasch described the 
superrotation in the viscous atmosphere of Venus from the planetary surface 
upwards. The zonal wind field thereby rigorously satisfied slip-free and 
stress-free boundary conditions. Due to the neglect of the viscous stress in 
the meridional momentum equation, however, the model is not entirely self 
consistent in that the meridional winds could not be made to vanish at the 
surface. Our results, which do not suffer from this deficiency nevertheless 
confirm the Importance of the meridional circulation for maintaining 
superrotation. 

What remains then to be shown is how, precisely, the flow Is built up or 

V 

superrotation is produced and if and how the atmosphere receive: angular 
momentum in this process? In answering these questions we shall ..^>cuss in 
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the first part of this paper slnpla mechanical analogies >:hat c Lucldate the 
basic process* Xn the second part we shall present a spectral aodel which 
describes in the linear regime the thermally driven tonally symnetrlc 
atmospheric circulation and superrotation. The approach is taken there to 
study the hypothetical spin up of a planetary atmosphere when differential 
solar heating is suddenly applied and follow the time evolution until a steady 
state is reached. Solutions obtained over a wide range of parameters thus 
provide Insight into the processes that establish and maintain superrotation. 
We find that, consistent with the mechanical analogies, two aspects are 
Important: With high viscosity the angular momentum for superr >tation is 

primarily derived from the change in the moment of inertia due to mass 
redistribution by the meridional ctrculatloii. With reasonable Low vlsocisity, 
however, this mass redistribution and the resulting pressure buLld-up in the 
polar region at lower altitudes causes the atmosphere to subrotate temporarily 
during the spin-up process so that the atmosphere can pick up the angular 
momentum from the planet giving rise to superrotation. 


LOCAL DYNAMICS 

We consider simplified equations to describe the zonally symmetric com- 
ponent of horizontal momentum conservation 


2 2 , 

2(iJ slnOU + — tgO + hn~s’ “ ““ 
p r * '0g^2 rp dQ 


( 1 ) 


2».ln0V + “!' tgS+_!' 1“ r, i!".o 

p r ® rcosO 30 r ^ 2 


( 2 ) 
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with U, zonal velocity 

V| Mtidional velocity 
p, aass density 

H0 viscosity coefficients in the vertical and horizontal iSirectlon 
' p, kinetic pressure 
(1)^, rotation rate of the planet 
0» latitude 

Assuming |v] « |u|and with a given pressure field, p(r, 9), the solutions of 
(1) and (2) yield In the Halt of small but finite viscosity (Hayr et al«, 
1981 ) the equation for gradient flow 


U ■ -w rcosO . . 

P l“„l P 


. “p / 2 2 2I ctp~3p 

+ J (ji r cos 9 -2- 

U.i I * r» p 30 


It reduces to 


U - 




2rpa)pSln0 


(|U| « |(. r|) 


and 



(3) 


(4) 


(5) 


in the geostroplc ami cyclostrophic approminations respectively. 

Formally, the above solutions represent a local balance between inertial 
forces and pressure gradients. They represent the simplest and perhaps iaost 
powerful concept for describing' superrotation. On most planets in the solar 
system the rotation axis is perpendicular to the orbital plane, the equatorial 
region is heated preferentially and the pressure tends to increase toward the 
equator (|£- < 0} at higher altitudes. The propensity toward superrotation in 

09 
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planetary ataraepheres Is thua explicable In terms of Equations ( V) through 
(5)> The largest superrotation rates are likely to occur at alt tudes where 
the viscous constraint In the meridional force balance la not Innortant, aiid 
the above solutions are thus approximately valid. An outstandln; example Is 
the work by Leovy (1973) where the large superrotation rate on Vmus was 
attributed to cyclostrophlc balance. 

Using the same horizontal pressure gradient, " ”0.3, fjr different 

planets, we note that the concept o£ local dynamical equilibrium produces a 

wide range of values 53 (Venus), 0.20 (Earth) and 0.0014 (Jupiter) which 

P 

match the observed average rotation rates within a factor of two. The zonal 
velocity Is here expressed In the form of rAo). 

AHGULAR MOMENTUM BUDGET 

The concept of local dynamical equallbrlum does not provide a connection 
between suphrrotatloh and source mechanism. Moreover, It Is restricted to 
that part of the atmosphere where, formally, viscosity Is unimportant and thus 
It cannot provide a mechanistic connection between the superrotatlng 
atmosphere and the rotating planet. 

The concept of angular momentum conservation Is fundamental for our 
understanding of super rotation. With the primary emphasis on elucidating the 
physics, we shall describe here a scenario of atmospheric superrotation which 
heurlstlcally evolves from fluid dynamical elements that can be readily 
understood to a more realistic configuration containing these elements. In 
support of this Interpretation we shall make use of mechanical analogies. 

As reference, a globally uniform atmosphere with mass density p (r) Is 

o 

chosen. It Is heated with a source, Qo(r), representing the global average of 
solar Input. Due to viscous Interaction this atmosphere corotates with the 
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angular velocity of the planet, Up. Without adding an external corque to the 
entire atnonphere’i'planet eyatem, we change the heat source to Q(r|6) ■ 

(^(r) + AQ(r,0) which In turn changes the density distribution to p(r,9) » 
PqCc) + Ap (r,0). The question we ash Is how AQ ealntalnK an a^aospherlc 
rotation rate w(r,0) " Wp + Aai(r,0) with a rigid shell componen: different 
froB that of the planet? 

In reality, the rotation rate of the planet also changes d>^\rlng the spin 
up process* For the purpose of discussing superrotation, howeve'^, we are 
justified to neglect this effect since the moment of Inertia of the atmosphere 
Is much smaller than that of the planet. 

Under conservation of angular momentum, L, the physical covdltions for 
uniform and nonuniform heating, representing end states of equilibria, are 
related 

\ (Q) » L^CQp) - AL, (6) 


where 


4L - yo) - I^(Q„) 


(7) 


represents the total amount of angular momentum that Is transferred from the 
atmosphere to the planet (positive) or visa versa (negative) during the spin 
up process. Subscripts a and p refer to atmosphere and planet respectively. 
Assuming Ap « p^, the explicit form of Eq. (6) is 




ApdiJ) - AL. 


( 8 ) 


8 



2 2 

With the definition di|» - r cos 0dt» the integrations ate taken tver the entire 
voluM of the atmosphere (with element dt). Equation (8) states that three 
components contribute to the conservation of the total angular iromentum: (a) 

atmospheric superrotation, (b) changes In the atmospheric mcnient of 
Inertia which are associated with changes In mass density, &p, and (c) the 
above defined exchange of angular momentum between atmosphere ai d planet, dL, 
during the transition from one (Qq) to the other (Q) equlllbrluir state. 

A simple thought experiment illustrates the nature of AL. starting out 
with an atmosphere characterized by Q^, p^, and o>^, corotating with the 
planet, we lower the heat source uniformly over the globe Q(r) ^ - 

AQ(r). Under mass conservation the atmosphere contracts, p(r) “ ” 

Ap(r) and its moment of Inertia decreases. During this change, the 
atmosphere tends to maintain its angular momentum and Initially there Is 
superrotation. This motion produces viscous Interaction with transfer of 
angular momentum to the planet, until eventually a new static equilibrium Is 
established In which atmosphere and planet corotate with the angular velocity, 
(i>. Of course, for this particular case the solution Is simply 


AL ■ Ip (o> - 0)^) > 0 


(9) 


with 


otlp+ foo'i*! 

TTi: y pdVl • 


where Ip refers to the planet's moment of Inertia. Considering Am - 0 another 
form of the solution can be obtained from Eq. 8, 
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(U) 


AL « 

This means that the atmosphere transfers niaount of angular rdiiantuB) 
equivalent to the entire change In Its ■onent of inertia; thus none reswlns in 
the atmosphere to aalntaln supcrrotatlon. 

PIROUETTE 

For a realistic atmosphere %»lth globally nonunlform sourci the confterva- 
tlon of angular mosMntum Is complicated, as will be shovRi in mnierlcal simula- 
tions of the spin up process* To provide a fundamental understanding of 
superrotation it is therefore hejpful to discuss first a simple model which 
proves to be valid for very large kinematic viscosities* 

The process is Illustrated in Figure 1. On a planet when the equatorial 
region is heated preferentially the compressible viscous atmosphere expands. 
Thus large scale motions in the form of Hadley cells are set u i (e.g* Lorenz, 
1967, Leovy, 1973, Gierasch, 1975) that effectively transport mass from the 
equator toward the poles closer to the rotation axis* As a result the 
atmospheric moment of inertia decreases 




AQdij) < 0 


( 12 ) 


which tends to decrease the angular n»raentum* Under Newton's law such a 
change is met with inertia and the angular velocity must Incre ise (super- 
rotation) to conserve angular momentum. We assume that during the initial 
phase of the spin up process the atmosphere superrotates all tne way down to 
the planetary surface and "reacts” by transferring, through a viscous torque, 
angular momentum to the planet. 

Unlike the case for a uniform source where this momentum transfer cor- 
responds to the entire change in moment of inertia (Eq* 11), only part of it 
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li tr«ntferr«d to the planet (&L>0) under nonuntform heating; this super- 
rotation can he sui6t'j*ined , which Is restated in oathematlcal fora 

<Au>I^ = d<|/ “ -Wp J j^pdiJ) -AL > 0. (13) 

Apart from the poleward transport of mass (term containing Ap In Eq. 13) 
which is driving iuperrotation* the meridional winds are also Irporf-sn'-. in 
•'Isolating'' the super rotating atmosphere from the planet. The < ssociated 
Inertial forces (which change sign as one vertically transverset a Hadley 
cell) are In balance with the viscous forces from the zonal wlnt' field, so 
that In equilibrium superrotatlort can be maintained with stress free 
boundaries. Solar heating and the resulting meridional clrcula! Ion thereby 
continually supply the energy that Is expended through viscous dissipation. 
During the spin up process (non-equlllbrlum) the Inertial force-, effectively 
limit the transfer of angular momentum to the planet. 

The magnitude of this momentum transfer is proportional to the viscosity 
and It has the effect of damping superrotation. A more viscous atmosphere Is 
tighter coupled to the planet which tends to constrain it into corotation. 

In support of this Interpretation a mechanical model is presented which 
Is closely related to the spin up In a pirouette. This analogy is crude by 
comparison with the dish pan experiments (e.g. Schubert and Whitehead, 1969) 
providing a fluid dynamic simulation of superrotation. However, it contains 
some of the Important elements to describe the process and has the advantage 
that the physics Is more transparent. 

The discrete phases of the pirouette experiment are illustrated in Figure 
2. Sitting on a rotating turntable, representing the solid planet, is a 
mechanical device referred to as "figure" which represents the atmosphere; its 
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configuration and aoaent of inartla (tilth aasaas In ahadad areaiS) can be 
changed by pulling a etrlng. The angular velocity of the figure ia achenat-* 
Ically ahown as a function of tiee. tn the flrat picture fraee the atring is 
relaxed and no external force ia applledn the figure haa the naxietli moment 
of inertia attainable and corotates Kith the turntable, which is the analogy 
of an atmoaphere under globally uniform heating* 

In the second frame the atring ia pulled. Expending energy againat 
gravity, Inertial forcea (e*g« centrifugal and Coriolis) and friction, the 
figure atraightens out, ita moment of inertia decreaaes and it superrotates. 
Locked in this position, the new figure would continue to superrotate indefin- 
itely if there were no frictional Interaction or momentum excharge with the 
turntable (thin solid line in the velocity diagram). Under realistic 
conditions, however, friction causes a slowdown (heavy solid lire in the 
velocity diagram) which eventually would lead to corotatlon. Tc prevent this 
the figure must be spun up periodically. This Is accomplished by relaxing the 
string so that the figure assumes its original moment of Inertia. All of the 
potential energy and p xrt of the kinetic (Inertial) energy previously expended 
are thereby released. In that position the figure subrotates and picks up 
angular momentum from the turntable. By pulling the string the figure is then 
spun up again, and a new cycle begins. 

The rotation rates in Figure 2 illustrate the steady state condition 
during an experiment cycle and Indicate that the net result of the procedure 
Is superrotation. A simple analysis of the mechanism can readily verify 
this. 

The equation of motion has the form 

(las) * - ^ las, 

la 
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d 

dt 


( 14 ) 


ORIGINAL PAGE IS 
OF POOR QUALITY 


where T(t) is the time varying moment of inertia, m is the difference between 
the angular velocities of the figure (with truss m) and the turntible, and a is 
the friction coefficient. Over the experiment cycle, in steady ^tate, the net 
amount of momentum exchange with the turntable must he zero 



d 

dt 


(lo)) dt = - 




(15) 


and a >2 being the average angular velocities during the two phases of the 
cycle and tj, and t 2 » I 2 being the corre'tponding time periods 
(t » tj^ -f and moments of inertia. Equation (15) becomes 


4 “1 




“2 "2 


= 0 . 


(16) 


It follows then immediately that the average superrotation rate is 



I 

X 


m. 



— SL 

T a 



- MJ2) O - -J-) > 


0 


(17) 


which is positive, since > 0 and < l 2 » The an,^ular vel- 

ocities and t0j^2 are defined in Figure 2. The figure superrotates because 
of the cyclic change in the moment of inertia. Energy is thereby continually 
fed into the system where it is dissipated through friction. The super- 
rotation rate increases inversely proportional to friction. 

Relative to the initial state of corotation (with “q)> the subsequent 
state of superrotation has a lower (average) angular momentum. The differ- 


13 



ence, Ai > 0 (Eq. 13) Is transferred to th'e turntable during thi spin 
process. It In'creases with friction attH decreases superrbratlot. . 

In sumitary, we note the a'dilogles with the superro'tatli^ at mosphiere. the 
oscillatory niotlons within one cycle of the plrottl^tte experiment simulates the 
Hadley cell. The higher the frequency Is In the experiment eye e, the shorter 
Is the turnover time in the Hadley cell and the more energy Is oltig through 
the superrotation engine. In the former, mass and momentum excl angd take 
place sequentially, while in t^te latter they occur Simultaneous y in the upper 
and lower legs of the circulatloh. The external force on the st ring repre- 
sents the horizontal pressure gradient due to differential heat ng. These 
forces drive the radial motions of the pirouette experiment and the meridional 
motions In the atmosphere, thereby controlling the moments of 1 ertla in both 
systems. They are primarily balances by Coriolis and centrifug..! forces which 
provide the link with friction and viscous dissipation where ehi*rgy is 
ultimately dissipated. Thus closing the energy chain we obtain a basic under- 
standing of Superrotai Ion which is consistent with the interpre atlons of 
Leovy, 1973 and Giera.' ch, 1975. 

One easily sees that the concept of this mechanical analog ' can als6 be 
applied to Individual particles of air moving with a Hadley cel . In a fluid, 
at any given time, every phase of this motion is equally repres -hted and thus 
justifies the time averaging procedure aS a means of representl ig the dynamic 
system in steady state. 

DOUBLE PIROUETTE 

The problem with the above analogy is that the entire atmosphere is 
portrayed as a body which can be deformed horizontally but not /ertlcally. 

Such a picture is essentially valid when large viscosity provld->s enough 
rigidity. Under more realistic conditions with reasonably low viscosity. 
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however, the situation Is note complicated. 

A logical extension of the diechanlcal model Is the "double pirouette" 
experiment which is significantly different from the earlier slrgle pirouette 
experiment in that Is has an additional degree of freedom for the vertical 
momentum distribution. Instead of one mechnical device, two art- used on top 
of each other, separated by a friction plate (Figure 3). Pulling the string 
at the upper figure (F^) supplies the energy. As Is stretched and Its 

moment of inertia decreases, the lower figure (Fj^) is compressed and Its 
moment of inertia Increases. Springs in Fj^ may provide a restoring force due 
to the compressibility of the atmosphere. 

The analytical solution for this thought experiment is straightforward. 
But depending on the conditions, the results may differ significantly, and a 
detailed discussion would go beyond the scope of this paper. Adopting here a 
simplified description, we consider the corotating (tu^) double pirouette with 
moments of Inertia and representing the atmosphere under globally 

uniform heating. The string is pulled during the "energisation phase" and the 
moments of Inertia change to I®, it Is relaxed during the "dissipation 
phase" and the moments of inertia return back, to I^, I^. Subscripts u and 1 
refer to the figures and Fj^, while superscripts e and d Indicate the 
energization and dissipation phases respectively. During the energization 
phase, the time average rotation rates (relative to o)^) are > 0 and 
w® < 0 in F^ and Fj^ respectively. 

Conservation of angular momentum with a system of two coupled equations 
corresponding to (15) yields 


d 
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where 

+ a 

Y 

a^, friction coefficient between turntable and pirouette device 
a, friction coefficient between both pirouette devices. 

Since I® < (5® > f) and I® > if (S)? < 0), it follows from (18) and (19) 

that the average rotation rates for both figures are positive 

“u • T ^^u *^*' ■*■ ^ ° (20) 

^1 “ T (“l ^ t**) > 0 (21) 

Examlng the angi lar momentum budget (Eqs. 8 and 13) 

0 < wl = S 1*^ + 3 , if - oj (I® - + I? - if) - &L (22) 

o uu 11 o u u 1 1' ' ' 

6 d 6 d 

produces an interesting result. Suppose we chose I^ ■ 1^^ and 1^^ * I^ which is 

the case Illustrated in Figure 3. Under this condition the total change in 

6 d 6 d 

the moment of inert! i is AI = (I^ - I^ + I^ - 1^^) ■ 0. Knowing from (20) and 
(21) that there is superrotation, Eq. 22 then implies that the double 
pirouette device is capable of picking up angular momentum, AL < 0, from the 
turntable during the transition from (uniform) corotation to super rotation. 


16 


This result Is completely different from that obtained with the single 
pirouette experiment where the decrease in moment of inertia caiises 
superrotation and, additionally, permits angular momentum to be transferred to 
the turntable during the spin up process. Depending on the amount of 

friction, the total change in moment of inertia, AI, and or the (initial) 

exchange of • angular momentum with the turntable (planet), AL, m.iy thus con- 
tribute to the angular momentum for superrotation. Under atmos;>heric 
conditions with reasonably low viscosity it is the latter which is usually 
more Important. 

The particular form that was used here to energize the double pirouette 

experiment is not unique. An external force applied to the low<*r figure would 

lead to virtually the same conclusions. 

On the basis of this analysis which is consistent with the results from 
our numerical simulation of atmospheric superrotation in the second part of 
this paper, the follov.'ing conclusions are drawn: 

Viscosity. In the case of identically zero viscosity or friction, the 
equations (20) and (21) cannot be applied. But assuming that initially the 
mechanical device (or the atmosphere under uniform heating), somehow, 
uniformly corotates, t ne can easily show that (u^ > 0 but < 0. The angular 
momenta in regions of super-and sub-rotation add up to zero implying that the 
average change in the rotation rate is also zero. Superrotation cannot be 
understood without vlrcosity. To sustain this motion the system must be 
capable of converting kinetic energy back Into thermal energy. 

Mass Transport. Without an external torque, conservation of angular 
momentum requires that the moment of Inertia must change in order to sustain 
superrotation. This is accomplished through mass redistribution by the 
meridional circulation of the atmosphere Induced by solar heating and pressure 
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gradients, or by the oscillatory motions due to extermal forcing in the 
prlouette experiments. To understand superrotation it Is Important that the 
change in moment of inertia is not maintained statically but through transport 
processes. The mass carrying meridional circulation of the atmosphere and the 
oscillatory motions of the mechanical analogies continually supply the 
energies that are expended by viscous dissipation associated with super- 
rotation. 

Momentum exchange with the planet. During the transition from corotation 
under globally uniform heating to superrotatlon under differential heating the 
atmosphere exchanges angular momentum with the planet. In simplest form, one 
expects a behavior similar to that of the single pirouette experiment where, 
through viscous interaction, angular momentum is transferred to the turntable 
thus slowing down superrotatlon. However, under more realistic conditions 
with reasonably low viscosity, a secondary circulation develops temporarily at 
lower altitudes and tlie analogy with the double pirouette experiment is 
appropriate. Initially, during the spin up process of the atmosphere or 
during the energization phase of the mechanical analogy, subrotation prevails 
near the surface, so that the systems can actually receive angular momentum 
from the planet or the turntable. 

Analogy with subrotation. The double pirouette experiment represents a 
mechanical model for nn atmosphere, preferentially heated at the equator. Our 
simplified description of this analogy, aimed toward a basic understanding of 
superrotatlon, however, cannot be adopted for subrotation. In reality, the 
changes in the moment of Inertia due to the oscillatory motions of the pir- 
ouette device are continuous and the Coriolis force must be considered. Under 
this condition the analysis is more complicated but shows that external com- 
pression against the action of the Coriolis force indeed leads to sub- 
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rotation* This analogy was used to Illustrate the subrotation on the sun 
whose atmosphere Is preferentially heated at the poles (Schatten et al., 

1981). Analogous to the fluid (e*g. Eq. 3), the mechanical experiment Is not 
symmetrical for large values of super- and sub-rotation when the centrifugal 
force becomes Important. 

FLUID DYNAMIC MODEL 

In the framework of a spectral model, we shall concentrate here on the 
lowest order harmonic or rigid shell component of superrotation which Is 
unique In that It alone pocesses a global angular momentum; the angular 
momenta of the higher order vector spherical harmonics are identically zero. 

A subsequent paper will discuss the problem of differential rotation and the 
acceleration of the equatorial region (Hide, 1971) by extending the analysis 
to higher order modes under consideration of mode coupling and viscous flow 
(Mayr and Harris, 1981 ). 

In the following sections the model Is presented, and boundary conditions 
employed at the planetary surface are discussed. Steady state solutions are 
then considered followed by a presentation and discussion of time evolution 
simulations. The thrust of this study Is toward a general understanding of 
superrotation rather than to provide quantitative results for a specific 
planetary atmosphere. 

To describe atmospheric superrotation, a three dimensional spectral model 
(Mayr and Volland, 1972; Harris and Mayr, 1975; Mayr et al., 1978, referred to 
as MH) Is adopted whose properties are summarized here: (1) Energy advectlon 

and convection associated with large scale circulation are considered. (2) 
Heat conduction and dissipative processes due to Ion drag and viscosity In 
both the zonal and meridional force balances are Included. (3) Vertical and 
horizontal mass transport Including momentum transfer between the major 


19 



species of the atmosphere are considered In self consistent form. (4) 
Perturbation theory Is applied relative to a globally uniform atmosphere which 
corotates with the underlying planet. The physical quantities Including 
temperature » densities and winds are represented In an expansion of spherical 
harmonics and time functions such as Fourier harmonics, which effectively 
separate latitude and time dependences from altitude variations. (5) The 
equations of energy, mass and momentum conservation are Integrated from the 
planetary surface up to the exosphere. At the lower boundary, the vertical, 
meridional and zonal velocities relative to the rotating planet are set to 
zero, and adiabatic conditions are assumed to prevail. For the upper boundary 
conditions, the properties of molecular heat conduction and viscosity require 
that the vertical gradients of temperature and horizontal winds vanish. 

In the zonally synmetrlc circulation (m*0) we consider the "fundamental 
mode" of superrotation 



- (r) P° (0) - X 2 (r) V 2 (3 sln^O - 1) 



- V, (r) ( 0 ) = V AwrcosG 

111 <P 


(23) 


(24) 

(25) 


where AQ 2 Is the external heat source, and AT 2 , Ap 2 » W 2 , V 2 , are the 

resulting perturbations in temperature, density and velocity components in the 

vertical, meridional and zonal directions respectively. Standard notations 

0 +0 ^0 

are used for the scalar- (P 2 ) and vector- ( 82 , C^) spherical harmonics. Due 
to the orthogonality of spherical harmonics, the contributions to the global 
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angular nomantum budget (Eq. 8) from the zonal wind field and thi't mass 

redletrlbutlon (or changes in the moment of Inertia) are 

/ p (r) r V C^r cos 0 dT s/p X V C° C? dt 
o ^ n n J n n 1 


■ (Cj dr «yp^Au>rco80C°dT 


(26) 


Apart from the corotating globally uniform atmosphere, 
mode or the rigid shell component of rotation (angular 
latitude) has non-vanishing angular momenta due to Au and Ap 2 > For higher 
order modes which represent differential rotation (variable angular velocity 
as a function of latitude) the global angular momenta are identically zero at 
each height. 

We emphasize again that in this paper no attempt is made to simulate a 
particular scenario of atmospheric (rigid shell) superrotation which would 
require a thorough analysis of the energy sources (due to e.g. solar 
radiation, auroral particles and Joule heating) and external forces (due to 
e.g. electric fields of solar wlnd-magnetospherlc origin) that are 
Important. For illustrative purposes we chose a global average atmosphere 
with Tq, which corotates with the planetary rotation rate representing 
the conditions on Venus. To drive superrotation, a specific heat 
source AQ 2 is adopted which peaks at an altitude of 50 km. With these Inputs 


dx 

(27) 

only the fundamental 
velocity Independent of 


21 


fixed, the height Independent vertical eddy diffusion coeffecient K and the 
corresponding coefficient of eddy viacoaity are varied to examine the effects 
on atmospheric superrotation* Since we consider here only the fundamental 
mode or the rigid shell component of superrotation the horizontal eddy dif- 
fusion coefficient does not enter the analysis* 

BOUNPARY CONDITIONS 

Disregarding the complexities due to composition, eight boundary 
conditions are required to solve the ordinary differential equations that 
describe the vertical structure of temperature, density and velocity fields* 

Of these, seven conditions are firmly established: (a) the vanishing vertical 

and horizontal (zonal and meridional) velocities at the lower boundary, (b) 
the requirement that the energy equation is satisfied under neglect of heat 
conduction at the lower boundary and (c) the assumption that the heat ' 
conduction-and momentum-fluxes vanish at the upper boundary* 

In steady state, one must further require for the rigid shell component 
of superrotation that the viscous torque or the momentum transfer vanishes at 
the planetary surface r^ 

o 

In Implementing the stress free boundary condition (28) it is important to 
realize that mass redistribution (and the resulting change in moment of 
Intertla) is one of the three elements that balance the steady state angular 
momentum budget (see Eq* 8)* This mass redistribution also produces a 
pressure variation at the surface, Ap 2 > which, a priori, is not known* Also 
contributing to the momentum budget is the atmospheric momentum exchange with 
the planet, AL, during the spin up process which is related to the surface 
stress. We require, therefore, that the surface ptassuri* or density 
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aaplltudes asiociated with P2** Are chosen such that the stress free boundary 
condition (Eq* 28) Is satisfied in steady state 

(iPj)r -0 (29) 

0 o 

In our computer model an iterative procedure is used combined with binary 
search to enforce condition (29)* 

It turns out that for a slowly rotating planet such as Venus this 
rigorous approach to the boundary value problem produces results reasonably 
close to those obtained with the (redundant) assumption of hydrostatic 
equilibrium^ This indicates that without constraint the system tends to find 
its natural state and In turn justifies using the hydrostatic condition for 
the time dependent solution. 

The concept of Influencing the angular momentum budget by varying the 
surface pressure distribution proves valuable In still another way. Thus a 
computer experiment can be performed In which the surface pressure 
ulstrlbutidn Is chosen so that the atmospheric angular momentum Is Invariant 
during the transition from corotation to superrotation 

(Ap2)r AL - 0. (30) 

o 

STEADY STATE 

For a fixed heat source, A Q^, which peaks at the equator, Figure 4 shows 

g 

computer results of rigid shell superrotation In steady state ( ~ ^) • 3Vo 

different values of the vertical kinematic viscosity corresponding to eddy 
diffusion coefficients of K ■ 10^ and 10^ are chosen. Solid lines represent 
realistic conditions for which the viscous stress or the vertical gradient of 
the zonal velocity is forced to zero at the planetary surface (Eq. 29). Under 
this condition, rigid shell superrotation prevails for both viscosities 
throughout the atmosphere. 
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In analyzing the angular BOBentun budget we evaluate the Integrals of 
^uperrotatlon and Bass redistribution f roe Eq. 8^ or Eqs* 26 and 27 

S •‘J’tui(r)p^v^coB^QiT (31) 

H - -Wg^p2(r) P2(0)r^cos‘OdT (32) 

M 

where AL •• M-S. The ratios 8 " are given in Figure A and reveal that In the 
case of large viscosity the atmosphere has lost angular aomentum to the 
planet AL > 0; while for the smaller value of viscosity the atmosphere has 
gained angular momentum from the planet AL < 0 during the spin up process. 

Retaining the Input parameters, we show for comparison a second set of 
solutions (dashed lines) which are artificial and diagnostic In nature. Here 
the surface pressure amplitude Is adjusted (Eq. 30) so that, under 
differential heating, AL«0. For large viscosity superrotation prevails all the 
way down to the surface. In that state the velocities are consistently larger 
than those under the stress free boundary condition (solid line), and a 
surface torque > 0 tends to slow down the excess superrotation while 

transferring angular I'lomentum to the planet. Obviously, such a condition 
cannot prevail In ste idy state. Initially, during the transition from 
corotation under uniform heating to superrotation under nonuniform heating the 
atmosphere tends to conserve Its angular momentum. Solutions for AL*0 are 
therefore Interpreted as heuristic representations of the transient state 
which, combined with the true steady state solutions “ 0, may serve to 

Illustrate the mechanics of the spin up process. 

For small viscosity the situation Is reversed. With AL ■ 0 the 
velocities are now consistently lower than those under the stress free 
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boundary condition. The atmosphere actually subrotates (minus) at lower 

altitudes. In this state a negative surface torque ^ ^ Induces the 

o 

atmosphere to pick up angular momentum from the planet which contributes to 
the larger superrotation. 

TIME EVOLUTION 


Steady state solutions contain no Information on causality and therefore 
cannot afford a complete understanding. An example Is the Hadley cell: Mass 

transport In; the upper leg of the circulation away from the source and in the 
lower leg of the circulation toward the source are In balance. But In the 
transient process of establishing the circulation, it Is precisely the 
Inbalance of mass flow which is building up the pressure in the polar region 
at lower altitudes, and It Is that pressure force which keeps the atmosphere 
In motion. 

We shall find that in principle the situation is similar for 
superrotation. There is, however, an important additional complication in 
that the atmosphere is coupled to a planet with virtually unlimited angular 
momentum. While the atmosphere's mass budget that governs the Hadley cell 
formation is closed, its angular momentum budget which governs supertotation 
Is not. 

In describing the spin up process we assume that up to the time t » 0 
the atmosphere is heated with a globally uniform source, AQ 2 ■ 0, and 
corotates with the planet, Au ■ 0. A globally nonuniform source is introduced 
as a step function 


AQ 2 (t) - 0; t < 0 
A (t) - const.; t > 0. 


(33) 
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Using L&place transformations, the time dependences of the heat source (33) 
and the physical variables X(t) (Eqs. 23, 24, 25) are described in the form 

m 

A (f) - / e"^^ A (t) dt - i (34) 

and 

q^ (f) x(f) - ^ X(t) dt, (35) 

respectively* With the time dependence e^^, f being a variable time constant, 
the MH model is used to compute the transfer function x(f) for unit heat 
Input. The Integral equation (*^5) is solved with the Gauss-Legendre technique 
of Bellman et al, 1966* Using six Gaussian points, the calculations are 
carried out over various timescales to provide the necessary resolution. 

For the same input values that produced the results of Figure 4, the time 
evolution of atmospheric superrotatlon is simulated and is given in Figures 5, 
6, and 7. For large viscosity the atmosphere superrotates at all altitudes 
throughout the spin up process. Immediately above the surface (2 km altitude) 
a transient velocity inaximum developes. In time, this velocity returns to a 

a 

small value near zero. The corresponding (transient) surface stress has the 
effect of decelerating and accelerating the atmosphere and the underlying 
planet respectively. That momentum transfer is also reflected in the steady 
state superrotatlon rates throughout the atmosphere which are lower than those 
during^ the initial phase of the spin up process. This confirms our earlier 
interpretation of the solutions in Figure 4. Judging from the time evolution 
of the height integrated superrotation rate we find that the atmosphere spins 
up like a solid body and transfers angular momentum to the planet, in analogy 
with the single pirouette experiment (Figure 2). 

For the smaller viscosity the results are entirely different (Figure 
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6). The Increased tendency tpyard geostrophy, combined with the Increase of 
surface pressure toward the poles (due to nerldlonal oass transport) Induce 
the staosphere to subrotate temporarily during the initial phase of the spin 
up below a critical level at about 40 km< Near the surface, the resulting 
viscous shear thus allows the atmosphere to pick up angular momentum from the 
planet which Is transported upward and, In time, during spin up, contributes 
to the large superrotation rates that eventually prevail at high altitudes In 
steady state. Leovy (1973) has suggested this process by Invoking the concept 
of geostrophlc adjustment. The phase of this wavelike phenomenon propagates 
downwards while energy and momentum are propagating upwards. We believe that 
below the critical level the angular momentum is carried upward by viscosity 
and gravity waves. Above that level, however, the downward propagating mode 
of viscosity waves may be more important while gravity waves continue 
traveling upwards until the spin up process is completed. The location of the 
critical level depends on the source distribution, and moves downwards as the 
viscosity Increases. 

The time evolution in the height distribution of mass flow is shown in 
Figure 7. Unlike the case with large viscosity where the change in moment of 
Inertia (Eq. 32) is important for the angular momentum budget (U 1), the 
rotational component (Eq. 31) dominates for R>1. Thus, a short period (2 
months) after energization the inertia of the system still tends to conserve 
angular momentum, and the regions of super- and sub-rotation are almost in 
balance. This condition is heuriscically described with Che AL » 0 solution 
In Figure 4b. As time progresses the atmosphere accumulates, through viscous 
stress, an Increasingly larger amount of angular momentum from the planet 
until eventually. In steady state, a stress free boundary condition is 
established and superrotation, Aa>>0, prevails all the way down to the 
surface. 

Considering that for low viscosity a large part of the atmosphere's 
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rotational energy comes from the planet, a scale analysis yields 

T - A 0) 

9 Ao) 

where the symbols K, and t stand for the viscosity, the vertical velocity 
gradient near the surface during the spin up process and the spin up time 
respectively. A decrease In viscosity is overcompensated by an Increase In 
the gradient and spin up time so that the atmosphere can receive an 
Increasingly larger amount of angular momentum. The meridional circulation 
transports angular momentum upwards and Its time response to the solar 
radiative Input Is almost Instantaneous. By comparison, the viscous flow of 
angular momentum from the planet is slow and slows further down as the 
viscosity decreases. Thus with decreasing viscosity, larger velocity shears 
must !» built up and a longer spin up time is required before the viscous flow 
can match the convective momentum transport. 

The dependence of atmospheric superrotation on viscosity is illustrated 
In more detail with Figure 8. For a fixed heat source and a wide range of 
eddy diffusion (viscosity) coefficients, the velocities are shown at various 
altitudes. To characierlze the angular momentum budget, the ratio between 
mass and rotational momenta, R, is presented. We identify two regimes: (a) 

For lower viscosities and R < 1 (AL < 0) the atmosphere receives an 
increasingly larger amount of angular momentum as K decreases. The velocity 
Is inversely proportional to viscosity, (b) At the higher viscosities K > 2 x 
10® the atmosphere is losing angular momentum to the planet and the velocity 
decrease Is more precipitous. For AL > 0 and AL < 0 the analogies with the 
single and double pirouette experiments are appropriate, respectively. 

CONCLUSION 

The spin up process leading to superrotation is discussed for the linear 
regime of the thermally driven zonally symmetric circulation. Depending on 
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the st:«te of turbulenc>^. the atnosphere can go through two' different modes. 

For large viscosity (K > 2x10^, in this study) it behaves almost like a solid 
body, tightly coupled to the planet. The superrotation rates are small and 
the spin up times are short. The angular momentum of superrotation is 
"driven'* by the decrease in the moment of inertia due to mass redistribution 
by the meridional circtilation. During the spin up time the atmosphere 
transfers angular momentum to the planet. This process has an analogy in the 
single pirouette experiment. 

For lower values of viscosity there is an Important additional degree of 
freedom in the vertical momentum distribution, and the analogy with the double 
pirouette experiment is appropriate. Mass can more freely return from the 
poles in the lower leg of the meridional circulation. The change in the 
moment of inertia is therefore not as large. This return flow combines witli 
the increased tendency toward geostrophy to cause subrotation in the lower 
region of the atmosphere during the Initial phase of the spin up process. 
Through viscous shear the planet applies a torque to the atmosphere which 
propagates up and in time supplies the angular momentum for large super- 
rotation, consistent with a suggestion by Leovy (1973). The spin up times in 
this process are long. Below a critical altitude level which is about 40km 
for K“10^ the mOmentuni flow during spin up is apparently carried by upward 
propagating gravity and viscosity waves. For a critical viscosity of about 2 
X 10^, in this analysis, the atmosphere does not exchange angular momentum 
with the planet. 

It is again emphasized that no attempt was made here to describe 
differential rotation and equatorial Jets, such as observed on Jupiter and 
Saturn. By expanding the spectral model to include higher order modes, the 
concept of mode coupling is used to describe differential rotation in an 
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atmosphere with zonal and nerldlonal viscous stresses (>iayt and Harris, 

1981). Momentum Is thereby allowed to cascade from lower order modes (which 
"absorb” most of the solar energy) to higher order modes, and the feedback on 
the rigid shell component of superrotation Is considered. At altitudes above 
the center of the meridional circulation, where the specific heat Input Is 
relatively large, angular momentum Is effectively advected poleward, except 
near the equator where the meridional winds and the Coriolis force vanish and 
where the horizontal momentum transport must be carried by horizontal diffu- 
sion (viscous flow). This causes the atmospheric rotation rate to decrease at 
low latitudes and Increase at high latitudes, relative to the rigid shell 
component of superrotation; Inducing a Cencency for the equatorial region to 
subrotate. However, below the center of the meridional circulation and below 
the nBxlmum energy input, energy and angular momentum are advected from high 
to low latitudes. Thus the equatorial region can be accelerated to the extent 
that positive jets develop. Again, near the equator, angular momentum must be 
carried by horizontal diffusion; the equatorial minima In the angular 
velocities of Jupiter and Saturn being evidence for down gradient momentum 
flow toward the equator. 

In this model the tropospheres of Jupiter and Saturn are assumed to be 
convectively unstable with a superadlabatic temperature lapse rate accounting 
for global scale energy diffusion toward higher altitudes where radiation to 
space becomes important. Under this condition, the upward motions In the 
direct and indirect atmospheric circulations Imposed by insolation can then 
supply dynamic energy for "localized" equatorial heating which in turn drives 
equatorial superrotation. Preferential solar heating at low latitudes Induces 
the ordered structure In the atmospheric motions. Tlie efficiency of this 
process is determined by the (in)8tabillty of the troposphere, the energy 
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conversion In the Interior and the radiative loss from the upper troposphere 
and stratosphere. Assuming zonal symmetry and thermal forcing by solar 
differential heating, we believe that four factors are Important In 
understanding the atmosphere dynamics of these planets: (1) Thelc large size 
and rotation rate produce through mode coupling a multiple Ferrel-Thomson type 
circulation with upward motions at low latitudes. Due to the lower mass, size 
and rotation rate of Saturn Its equatorial Jet Is much broader tlian that of 
Jupiter. (2) Energy from the Interior Is diffusing radially outwards and 
maintains a supeadlabatlc lapse rate which provides the condition for trapping 
some energy that Is advected toward the equator In the lower leg of the 
meridional circulation. (3) Tlte eddy dlffusivlty Is sufficiently low to 
require large temperature and velocity gradients for the energy and momentum 
balances, the atmosphere of Saturn being perhaps less turbulent than that of 
Jupiter. (4) The observed temperature differences on both planets (Hanel et 
al. , 1981) cause the cloud layers, the tracers of wind systems, to be formed 
higher up below the tropopause on Jupiter where the upper and lower legs of 
the meridional cell produce a complex Interference pattern In the angular 
momentum budget but deeper down In the convection region of the troposphere on 
Saturn where meridional energy and momentum advectlon toward th equator 
dominate. 
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Figure Captions 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


I Planet with rotation axis penpendlcular to the orbital plane around 
the Sun. Schematic diagram Illustrating the energy deposition In 
the equatorial region, the Hadley cell and the resulting mass relo- 
cation toward the poles. Under the tendency to conserve angular 
momentum the dynamically Induced decrease In the moment of Inertia 
produces superrotation. 

: Single pirouette experiment as analogy for an atmosphere with high 

viscosity. A mechanical device (figure), Its moment If Inertia 
being cyclically changed by an external force (representing solar 
differential heating and pressure gradient). Is frlctlonally 
coupled to a rotating turntable. Averaged over an "experiment 
cycle" the figure superrotates receiving the angular momentum from 
the (average) decrease in moment of Inertia. During the spin up 
process, angular momentum, AL, Is frlctlonally transferred to the 
turntable which slows down superrotatlon. 

: Double pirouette experiment as an analogy for an atmosphere with 

(reasonably) low viscosity. Compared to Figure 2 a second figure 
provides the added degree of freedom for the vertical momentum 
distribution. By alternatlngly changing the moments of Inertia of 
both figures they superrotate In the time average. But this can 
now be accomplished without changing the total moment of Inertia, 
thus the device Is capable of picking up the angular momentum for 
superrotatlon from the turntable. 

: Computed vertical distribution for the rigid shell component of 

"superrotatlon", UcosG. Velocities are zero at the surface. In 
solid lines steady state solutions are presented with stress free 
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boundary conditions* In dashed lines an artlficlnl condition is 
shown for which the angular monenttn of the atmosphere Is 
conserved, AL - 0* It Is Interpreted as a heuristic representation 
of the transient state during s|[>ln up. The results are given for 
two viscosities corresponding to eddy diffusion coefficients of K “ 
10^ (Figure 4a) and K ■ 10^ (Figure 4b). Plus and ndnus signs 
Indicate super- and sub-rotation respectively. 

Figure 5: Simulation of the time evolution leading to superrotation for K * 

10^. The latitude dependent component of solar heating which 
causes superrotation is turned on abruptly as indicated in the 
dashed line. Note that the atmosphere super-rotates throughout the 
spin up, with the positive velocity above the surface providing the 
viscous stress for momentum transfer to the planet during that 
process. 

Figure 6: Results for the same heat source as in Figure 5 except that K > 

10^. The atmosphere at lower altitudes subrotates during the 
Initial phase of the spin up. The negative velocity above the 
surface is picking up a torque form the faster rotation planet 
which provides most of the momentum for superrotation. Note that 
In both results (Figures 5 and 6) superrotation prevails all the 
way down to the planetary surface in steady state. With low 
viscosity, the spin up times are longer and the superrotation rates 
are larger compared to the corresponding values for large 
viscosity. 

Figure 7i The results from Figure 3 are shown for mass flow, P^U, at various 
(approximate) times during the spin up .process. Npte that 
Initially, about 2 months after energization, the angular momenta 
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In regions of super- and sub-rotation alaost balance each other, 
thus justifying our Interptetatlon of the AL ■ 0 solutions In 
Figure 4. 

Figure 8: Dependence of the steady state superrotation rate on viscosity* 

Also shown Is the ratio R between the total angular momenta due to 
changes In the atmospheric moment of Inertia and rotation rate. 

For R » 1 superrotation Is primarily driven by the mass 
redistribution (Figure 4a and 5) with angular momentum being 
transferred to the planet during spin up. For R « 1 supercotatlon 
la primarily driven by the angular momentum picked up from the 
planet during spin up (Figures 4b, 6 and 7). For R > 1 (AL > 0) an 
analogy with the single pirouette can be made, while 
for R ^ 1 (AL ;C 0) the analogy with the more corapllcated double 
pirouette must be Invoked. For the critical value 

O 

of K >■ 2 X 10 the atmosphere does not exchange angular momentum 
trlth the planet. 
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